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We investigate a possibility to measure parameters of the Littlest Higgs model with 
T-parity at the ILC with y/s = 500 GeV. The model predicts new gauge bosons (Ah, 
Zn, and Wh), among which the heavy photon (Ah) is a candidate for dark matter. 
Through Monte Carlo simulations of the process of e + e~~ — + AnZn, we show that the 
masses can be determined with accuracies of 16.2% for mj H and 4.3% for mz H - In 
additionally, the vacuum expectation value can be determined with accuracies of 4.3%. 

1 Introduction 

The Littlest Higgs Model with T-parity is one of the attractive candidates for physics beyond 
the Standard Model because it solves both the little hierarchy and dark matter problems 
simultaneously. One of the important predictions of the model is the existence of new heavy 
gauge bosons, where they acquire mass terms through the breaking of the global symmetry 
necessarily imposed on the model. The precise measurements of their masses allow us to 
determine the most important parameter of the model, namely the vacuum expectation 
value of the breaking. 

At the Large Hadron Collider (LHC), it is difficult to determine the properties of heavy 
gauge bosons, because they have no color charge [2J. On the other hand, the ILC will 
provide an ideal environment to measure the properties of heavy gauge bosons. Heavy 
gauge bosons are expected to be produced in a clean environment at the ILC, which enables 
us to determine their properties precisely. We study the sensitivity of the measurements to 
the masses of heavy gauge bosons and the vacuum expectation value at the ILC based on a 
realistic Monte Carlo simulation. 



2 Simulation framework 



In order to perform a numerical simulation, we need to choose a representative point in 
the parameter space of the Littlest Higgs model with T-parity. Firstly, the model parame- 
ters should satisfy the experimental results of the precise measurements on the electroweak 
parameters [3]. In addition, the cosmological observation of dark matter relics also gives 
important information. For that reason, we consider not only the electroweak precision mea- 
surements but also the WMAP observation [4] to choose a point in the parameter space. At 
selected representative point, the vacuum expectation value(/) is 580 GeV, and the masses 
of the heavy gauge boson are 81.9 GeV and 369 GeV for mA H and mz H , respectively. 
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Process 


Cross sec. [fb] 


# of events 


generated 


after all cuts 


A H Z K -> A H A H bb 


1.05 


525 


272 


vvh — > j^fofo 


34.0 


17,000 


3,359 


Z/i — ► i/i/6fo 


5.57 


2,785 


1,406 




496 


248,000 


264 


ZZ — ■+ 


25.5 


12,750 


178 


— > z/i/66 


44.3 


22,150 


167 


7Z — > 766 


1,200 


600,000 


45 



Table 1: Signal and backgrounds processes. 




Figure 1: 
process. 



Diagram for signal 



Because sum of mA H and mz H is less than 500 GeV, mA H 
and mz H can be generated through the process e + e~ — > AhZh 
at = 500 GeV. Feynman diagram for the signal process is 
shown in Fig. [TJ Zh decays into Anh with almost 100% 
branching fractions. 

We have used MadGraph [5] to generate signal events of the 
Little Higgs model, while Standard Model events have been 
generated by Physsim [BJ. Initial-state radiation and beam- 
strahlung have not been included in the event generations. 

We have ignored the finite crossing angle between the electron and positron beams. In both 
event generations, hclicity amplitudes were calculated using the HELAS library [7], which 
allows us to deal with the effect of gauge boson polarizations properly. Phase space integra- 
tion and the generation of parton 4- momenta have been performed by BASES/SPRING [8]. 
Parton showering and hadronization have been carried out by using PYTHIA6.4 [9], where 
final-state tau leptons are decayed by TAUOLA [10] in order to handle their polarizations 
correctly. 

The generated Monte Carlo events have been passed to a detector simulator called JS- 
FQuickSimulator, which implements the GLD geometry and other detector-performance 
related parameters [TTJ . 



3 Analysis 

The simulation has been performed at y/s = 500 GeV with an integrated luminosity of 500 
fb -1 . The heavy gauge bosons, Ah and Zh , are produced with the cross section of 1.9 
fb. Since Zh decays into Ah and the Higgs boson, the signature is a single Higgs boson 
in the final state, mainly 2 jets from h — * bb (with a 55% branching ratio). We, therefore, 
define Ah^h — > Ah Ah bb as our signal event. For background events, contribution from 
light quarks was not taken into account because such events can be rejected to negligible 
level after requiring the existence of two &-jets, assuming a 6-tagging efficiency of 80% for 
&-jets with 15% probability to misidentify a c-jet as a &-jet. This ^-tagging performance was 
estimated by the full simulation assuming a typical ILC detector. Signal and background 
processes considered in this analysis are summarized in Table [TJ 

The clusters in the calorimeters are combined to form a jet if the two clusters satisfy 
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Vij < Vent- Vij 



is defined as 



Vij = 



2E i E j (l - cos%) 

UK. ' 



(i) 



where is the angle between momenta of two clusters, Euj\ are their energies, and E v - ls is 
the total visible energy. All events are forced to have two jets by adjusting y cut . We have 
selected events with the reconstructed Higgs mass in a window of 100-140 GeV. In order 
to suppress the vvh — > vvbb background, the transverse momentum of the reconstructed 
Higgs bosons (px) is required to be above 80 GeV. This is because the Higgs bosons coming 
from the WW fusion process, which dominates the vvh ~ » vvbb background, have pt mostly 
below W mass. Finally, multiplying the efficiency of double 6-tagging (0.8 x 0.8 = 0.64), we 
are left with 272 signal and 5,419 background events as shown in Table[TJ which corresponds 
to a signal significance of 3.7 (= 272/V5419) standard deviations. The indication of the 
new physics signal can hence be obtained at y/s — 500 GeV. 

The Ah and Zh boson masses can be estimated from the edges of the distribution of 
the reconstructed Higgs boson energies. This is because the maximum and minimum Higgs 



boson energies (E n 



i) are written in terms of these masses, 



(2) 



where /?z h (7^h) is the /3(7) factor of the Zh boson in the laboratory frame, while E^ip^) is 
the energy (momentum) of the Higgs boson in the rest frame of the Zn boson. Note that 
E* h is given as (Mf H + Mf % - M\ h )/(2M Zh ). 
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Figure 2: (a) Energy distribution of the reconstructed Higgs bosons with remaining backgrounds 
after the selection cuts, (b) Energy distribution of the Higgs bosons after subtracting the back- 
grounds. The distribution is fitted by a line shape function determined with a high statistics signal 
sample, (c) Probability contour corresponding to 1- and 2-cr deviations from the best fit point in 
the An and Zn mass plane. The shaded area in shows the unphysical region of ttlah + rnzn > 500 
GeV. 

Figure Ufa) shows the energy distribution of the reconstructed Higgs bosons with re- 
maining backgrounds. We obtained the signal distribution after subtracting backgrounds as 
shown in Figl2]Jb). The endpoints, -Emax and E m i n , have been derived by fitting the distri- 
bution with a line shape function which was determined by a high statistics signal sample. 
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The fit resulted in m^ H and mz m being 83.2 ± 13.3 GeV and 366.0 ± 16.0 GeV, respectively, 
which should be compared to their true values: 81.85 GeV and 368.2 GeV. Figure^c) shows 
the probability contour for the masses of An and Zn . Since the masses of the heavy gauge 
bosons depend on the vacuum expectation value, / can be determined by measurement of 
the gauge boson masses. We obtained / = 576.0 ± 25.0 GeV. The input value of / is 580 
GeV in our simulation study. 

4 Summary 

At the ILC with y/s = 500 GeV, it is possible to produce An and Zn bosons with a signal 
significance of 3.7-sigma level. By using the energy distribution of the Higgs bosons from 
the Zn decays, the masses of these bosons can be determined with accuracies of 16.2% for 
mA H and 4.3% for toz h - From the results, it turns out that the vacuum expectation value 
can be determined with accuracy of 4.3%. The details on this analysis can be found in (T2"] . 
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